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Abstract 

We discuss the Schwarzschild solution in the Dvali-Gabadadze-Porrati (DGP) model. We 
obtain a perturbative expansion and find the explicit form of the lowest-order contribution. 
By keeping off-diagonal terms in the metric, we arrive at a perturbative expansion which is 
valid both far from and near the Schwarzschild radius. We calculate the lowest-order contri- 
bution explicitly and obtain the form of the metric both on the brane and in the bulk. As 
we approach the Schwarzschild radius, the perturbative expansion yields the standard four- 
dimensional Schwarzschild solution on the brane which is non-singular in the decoupling limit. 
This non-singular behavior is similar to the Vainshtein solution in massive gravity demonstrating 
the absence of the van Dam-Veltman-Zakharov (vDVZ) discontinuity in the DGP model. 
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Extra dimensions have been successful in explaining the weakness of gravity [1-5]. They 
modify Newton's Law by introducing Kaluza-Klein modes. In the case where the extra volume is 
infinite, light Kaluza-Klein modes dominate even at low energies [6-8] , leading to a modification 
of Newton's Law at astronomically large distances [9-17]. However, this may be avoided in the 
case of more than one extra dimension [7]. Problems arise when the brane is given a finite 
thickness in the transverse directions [18-20] obscuring the zero-thickness limit. 3 

The case of a four- dimensional brane living in a five-dimensional space of vanishing cos- 
mological constant (Dvali-Gabadadze-Porrati (DGP) model [6]) appears to be plagued by the 
van Dam-Veltman-Zakharov (vDVZ) discontinuity [22,23] and has recently attracted some at- 
tention [24-31]. This is similar to four-dimensional massive gravity where one does not obtain 
agreement with Einstein's General Theory of Relativity in the limit of vanishing mass of the 
graviton. Vainstein [32] provided a solution to the apparent discontinuity in the case of a point 
source by suggesting that the discrepancy arises from the linear approximation to the full field 
equations which has a limited range of validity. This second solution reduced to the Schwarzschild 
solution in the zero graviton mass limit demonstrating the absence of the vDVZ discontinuity in 
this spherically symmetric case. 4 

Applying a similar procedure to the DGP model is not straightforward, because the non- 
linear field equations are too complicated to solve even in the spherically symmetric case of a 
point source. It was subsequently argued in [27] that the discontinuity can be avoided by bending 
the brane through a coordinate transformation. This bending can be large even if the matter 
sources are weak. In refs. [28-31] solutions were found interpolating between the regimes far from 
and near the Schwarzschild radius by keeping second-order bulk derivatives in the perturbative 
expansion. It was thus shown that in the decoupling limit one recovers the standard four- 
dimensional Schwarzschild metric. However, the method of solution in [28-31] only yielded the 
form of the metric on the brane and provided no information on the bulk dependence of the 
metric. 

Here we discuss a perturbative solution to the DGP field equations in the case of a point 
source. By employing an ansatz for the metric containing off-diagonal contributions, we arrive 

3 See [21] for a consistent, finite-thickness brane model. 
4 See [33] for problems associated with Vainstein's approach. 
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at a lowest-order approximation to the field equations which is solvable. The solution is found 
explicitly throughout its domain of validity (both near and far from the Schwarzschild radius), 
including the bulk. 

The DGP model [6] describes a 3-brane on the boundary of a five-dimensional bulk space E. 
The action is 

S = M 3 [ d 4 xdyV^Gn {5) + M 2 f d 4 x\/^n {4) (1) 

where 7£( 5 ) (7^(4)) is the five-(four-) dimensional Ricci Scalar. The solution to the linearized 
equations bares a striking resemblence to the vDVZ solution of massive gravity and shares the 
apparent discontinuity in the decoupling limit (M — > 0). Porrati [27] argued that this solution 
is only valid in a limited domain, as was the vDVZ solution of massive gravity, and breaks down 
in the regime 

/ mM 2 \ 1/3 



r 



~ rc ' rc= {l8^J (2) 
when a static spherically-symmetric source of mass m is considered. 

We seek to obtain a solution to the field equations for a static point source with stress-energy 
tensor 

Tab = m5° A 5° B 5 3 (x)5(y) (3) 
which is valid throughout the region r m < r < l/m&, where 

2m 



m 8ttM 2 

is the Schwarzschild radius and 



(4) 



mb = Wr (5) 

is a crossover scale between four- dimensional and five-dimensional behavior. In the decoupling 
limit, mi — > and r is unbounded from above. 
We choose the ansatz for the metric 

ds 2 = -e 2B(r ' y) dt 2 + e 2C ^5 ij dx i dx j + 2^(r, y)dx i dy + e 2D ^ y) dy 2 (6) 

To arrive at a set of compatible field equations, we shall keep first-order contributions in the 
diagonal components B, C, D, and up to second-order terms in the off-diagonal field A. It is also 
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convenient to introduce the notation 

1=V0 , = V) 2 (7) 

r 

where prime denotes differentiation with respect to r, the distance from the point source. A dot 
will be used for bulk derivatives (with respect to y). 

First, let us discuss the lowest order contributions to the field equations in the bulk. The yy 
component reads 

2di&(B + 2C) + dA&^di&t) - 80^8° <P) = (8) 
The mixed components are 

di(B + 2C) - d^djditp = (9) 
The spatial brane worldvolume components are 

(didj - 5ijd k d k )(B + C + D - 0) - S^B + 2(7) 

+d k (d k <f ) d l d 3 <j ) ) - ^(dVW) - ISij (dki&tdt&ft - dk^dtd"^)) = (10) 
and finally, the tt component is 

d i d i {2C + d - 0) + 3C + \ [dkid^did^) - djid^dk&cj))] =0 (ii) 

In terms of the field \1/ (eq. (7)), the field equation (8) becomes linear, 

B + 2C + ^ = (12) 
Then the mixed components (9) may be written as 

d& + &><t>djd^ = (13) 

whose general solution is 

tf(r,y) = ^ + P(r) (14) 
We shall make use of gauge freedom and demand (as well as \1/) be independent of y, 

= ^ = (15) 

This is true to lowest order; higher-order corrections will introduce a non- vanishing (and ty). 



The remaining field equations (10) and (11) also become linear. They read, respectively, 

B + C + D + ±V = (16) 
V 2 {2C + D + ^) + 3(7 = (17) 

where we used (15). 

Eqs. (12), (15) and (16) yield 

C=-!(B + tt) , D = -\B (18) 

Then eq. (17) becomes 

B + V 2 B = (19) 
whose solution is easily obtained after Fourier-transforming the worldvolume coordinates 

B(p,y)=B(p,0)e~ p y (20) 

The bulk behavior of the other fields is then found from (18), 

C(p,y) = -lB(p,0)e-vy-^(p) , D(p, y) = -\B{p, 0)e~" (21) 

Having obtained the functional dependence of all fields on y in terms of data on the brane (y — 0), 
we now turn to solving the boundary field equations. 

On the boundary (y — 0), the spatial brane components yield 

2M 3 (f) + M 2 (B + C) = (22) 

whereas the tt component is 

QM 3 C + 2M 2 d i d i C + 2M 3 d i d i <f) = -md 3 (x) (23) 

The first term in eq. (23) may be dropped in the regime of interest, p ^> = 2M 3 /M 2 . 
Eliminating C by using (18), we obtain 

M 2 V 2 (B + - 2M 3 V 2 = m5 3 (x) (24) 

Solving for B on the boundary, we find 

vn 

5(r,O) = -^+m b 0- i -^ (25) 



Using this and (18) to eliminate B and C from the other boundary eq. (22), we deduce 

3 TTt 

Differentiating with respect to r and using (7), 

5^-;M ,+ dk = ° (27) 

which can easily be solved for <f>' yielding, 



4 \ V r 3 
where r c is given by eq. (2). We also have 



4>'=- A m h r[\- \l + r i) (28) 



r = I(0f = m + 9 2 r L_ 1 + r l\ (29) 
r Vr ; 8ixM 2 r 2 8 b \ \ r 3 J K J 

Differentiating eq.(25) with respect to r, we obtain the form of the field B on the brane, 



B ' ( ''-°^83V-^''( 1 -V 1 + 5) < 3 °> 

Summarizing, eqs. (18), (28) and (30) provide the form of the metric on the brane. This solution 
is valid everywhere on the brane (in the regime r m < r < l/m^). The solution in the bulk is 
given by eqs. (20) and (21) where the fourier transform B(p, 0) is deduced from (30). 

Next, we examine the near and far regimes, separated by the crossover distance r c (eq. (2)). 

In the far regime (r > r c ), we have from (28), 

and so (using eq. (7) or (29)), 

m , „/„6 



* = =t- { "I + o(r-°/r b ) ^ (32) 

1287rM 2 r \r 3 v c/ ; J v ; 

Notice that \& is of higher order. Therefore, eq. (18) implies to lowest order 

C(r,y) = D(r,y) = -±B(r,y) (33) 

Using (30), we obtain 

^^"sSf' + ^i (34) 
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and after fourier transforming, 

%,0) = -^ {1 + o(pVpI)} (35) 
where p c ~ l/r c . The y dependence of the fields to lowest order is given by 

~ ~ ~ 171 

C(p,y) = D(p,y) = -\B{p,y) = -^ 2 e~ py (36) 
Now taking the inverse fourier transforms, we finally obtain for r > r c , 



= 24^( 1 + ° (r ' /r3) ) (37) 

This solution corresponds to that of the standard perturbative expansion. Notice that in the 
decoupling limit (M 3 — > 0), <f> diverges. But this is a limit beyond our approximation above, 
because corrections become infinite (r c — > oo, from eq. (2)). We next obtain the solution in the 
near regime which corresponds to the four-dimensional Schwarzschild solution and is valid in the 
decoupling limit. 

In the near regime (r < r c ), we obtain from eqs. (28) and (29), respectively, 



mr r (/ , v ,„ vl „ m 



2nM* ( 1 + °((^) 3/2 )} > $ = "Mft (l + o((r/r c ) 3 / 2 )} (38) 
In this case \l/ contributes to lowest order. Using (30), we deduce on the brane 

The y-dependence to lowest order is given by 

cm = a^+^i 

b(M) = w 1 "" < 40) 
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where we used (20) and (21). At y — 0, we recover the Schwarzschild solution and therefore 
agreement with the standard Newtonian potential of massless gravity. All fields are now non- 
singular in the decoupling limit M — > 0. Fourier transforming, in the regime r m < r < r c , we 
obtain from eqs. (38) and (40), 

B<(r,y) = (l--tan-^)(l + o((r/r c ) 3 / 2 )) 

8ttM r I t r J 

C<(r,y) = -^p |l-itan-^)(l + o((r/r c ) 3 / 2 )) 
8ttM r I t r J 

D<(r,y) = — ^- {l--tan-^)(l + o((r/r c ) 3 / 2 )) 
16ttM r I t r J 

0<(r,y) = (l + o((r/r c ) 3 / 2 )) (41) 

To conclude, we derived a perturbative expansion for the Schwarzschild solution in the DGP 
model. By keeping terms in the field equations which were second-order in the off-diagonal 
metric components, we arrived at an explicit form of the solution both on the brane and in 
the bulk. On the brane, our solution interpolated between the near and far regimes separated 
by the distance scale r c (eq. (2)). At distances below the critical radius r c , the perturbative 
expansion yielded the four-dimensional Schwarzschild solution, demonstrating the absence of the 
van Dam-Veltman-Zakharov (vDVZ) discontinuity [22,23]. 
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